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The fully nonlinear governing equations for spin | quantum plasmas are presented. Starting from 
the Pauli equation, the relevant plasma equations are derived, and it is shown that nontrivial quan- 
tum spin couplings arise, enabling studies of the combined collective and spin dynamics. The linear 
response of the quantum plasma in an electron-ion system is obtained and analyzed. Applications 
of the theory to solid state and astrophysical systems as well as dusty plasmas are pointed out. 
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There is currently a great deal of interest in investi- 
gating collective plasma modes [l], 0, 0, 0, H, H, 0, [1] in 
quantum plasmas, as such plasmas could be of relevance 
in nano-scale electro-mechanical systems 0, [T3, El, m 
microplasmas and dense laser-plasma s 11211 , and laser in- 
teractions with atomic systems 13, LLJ]. For example, 
Refs. [J and @, i, d used quantum transport models 
in order to derive modified dispersion relations for Lang- 
muir and ion-acoustic waves, while Shukla & Stenflo [151 ] 
ionvestigated drift modes in nonuniform quantum mag- 
netoplasmas. Moreover, it is known that cold quantum 
plasmas can support new dust modes 0, [13 ■ In Rcf. 
[8( it was shown that electron quantum plasmas could 
support highly stable dark solitons and vortices. Fur- 
ther examples of quantum plasmas and the range of va- 
lidity of their descriptions has been discussed recently 
in Ref. [HI- The above studies of quantum plasmas 
have used models based on the Schrodinger description 
of the electron. It is expected that new and possible 
important effects could appear as further quantum ef- 
fects are incorporated in models describing the quantum 
plasma particles. The coupling of spin to classical mo- 
tion has attracted interest in the literature (see, e.g., 

1111,11,11 m m m m m s m m oni). Much 



work has been done concerning single particle spin effects 
in external field co nfig urations, such as intense laser fields 
12, H, 24, 25, 26|, 27 1, and the possible experimental sig- 



natures thereof. However, there have also been interest 
in excitations of collective modes in spin systems, such as 
spin waves, in a wide scientific community. For example, 
in Refs. [3, H HH hydrodynamical models including 
spin was presented, and further theory concerning spin, 
angular momentum, and the forces related to spin was 
discussed in Refs. and [1(3]. Moreover, spin waves 
in spinor Bose condensates has recently been discussed 
in, e.g., Ref. [3lj]. The treatment of charged particles 
and plasmas using quantum theory has received atten- 
tion in astrophysical settings, especially in strongly mag- 
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netized environments 32j, [33]. For example, effects of 
quantum field theory on the linear response of an elec- 
tron gas has been analyzed [131, results concerning the 
spin-dependence of cyclotron decay on strong magnetic 
fields has been presented [35| . and the propagation of 
quantum electrodynamical waves in strongly magnetized 
plasmas has been considered [36j . 

In this Letter we present for the first time the fully non- 
linear governing equations for spin \ quantum electron 
plasmas. Starting from the Pauli equation describing 
the nonrelativistic electron, we show that the electron- 
ion plasma equations are subject to spin related terms. 
These terms give rise to a multitude of collective effects 
of which some are investigated in detail. Applications of 
the governing equations are discussed, and it is shown 
that under certain circumstances the collective spin ef- 
fects can dominate the plasma dynamics. 

We will assume that the electron wave function can 
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be written in the product form 'F 
where N is the number of particle states. Thus we will 
here neglect the effects of entanglement and focus on the 
collective properties of the the quantum electron plasma. 
Then the non-relativistic evolution of spin ^ particles, as 



described b> 
(see, e.g. [37 



the two-component spinor &( a ) , is given by 
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where a numbers the particle states, m e is the parti- 
cle mass, A is the vector potential, e is the magnitude 
of the electron charge, [Ib = —eh/2m e c is the electron 
magnetic moment, <f> is the electrostatic potential, and 
<j = (o"i, 02, (T3 ) are the Pauli spin matrices, represented 
by 
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By introducing the decomposition of the spinors ac- 
cording to ^(a) = ^ n (a) exp(z5( a ) IK)<P( n ) , we may derive 
a set of coupled fluid equations [37( for the densities 
ni a ), the velocities vr a ) = (l/m e ) {VSr a ) — ihip^'Vip) + 



2 



(e/m e c)A, and the spin vectors S( a ) = (%/2)(p\ a -\tr<p^ a ) 
(where <p( a ) is the 2-spinor through which the spin ^ 
properties are mediated). 

Next we define the total particle density for the species 
with charge q according to n e = Y^,Ux)=i Pa n (a)> where 
p a is the probability related to the wave function ^( a )- 
Using the ensemble average (/) = J2 a Pa( n {a)/n)f for 
any tensorial quantity /, we define the total electron fluid 
velocity for charges V e = (vf a \) and the total electron 
spin density S = From these definitions we can 

define the microscopic velocity in the electron fluid rest 
frame according to itf( Q ) = V( Q ) — V ei satisfying (w^) — 
0, and the microscopic spin density <S( Q ) — s (a)~ S, such 
that (5 (q) ) = 0. 

We then obtain the conservation equations 



d t n e + V • (n e V e ) = 0, 
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mn e (d t + V e ■ V) V e = -en e (E + V e x B) 
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respectively. Here we have added the electron-ion colli- 
sions C e i, denoted the total quantum force density by 
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V-{n e Sym[(V^)®((V5 ( a Q) ))]}, (6) 



where Sym denotes the symmetric part of the tensor, and 
defined the nonlinear spin fluid contribution by 

Ms = -^-S x [d a (n e d a S)} + —S x [d a (n e (d a S (a) ))} 



rn, 
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where n e = TO e n e [(iO( a ) (g)W( a )) — l(w^ a j)/3] is the trace- 
free anisotropic pressure tensor (I is the unit tensor), 
P e = m e n e (w^ a yi is the isotropic scalar pressure, Z = 
(V5 a ) <g) (VS°) is the nonlinear spin correction to the 
classical momentum equation, I! = ((V5( Q ) a )(8)(V l S° Q j)) 
is a pressure like spin term (which may be decomposed 
into trace-free part and trace), K = n e (w;( Q ,)(g)<S( Q ,)) is the 
thermal-spin coupling, and [(V <g> B) ■ S] a = {d a B b )S b . 
Here the latin indices a, b, . . . = 1,2,3 denote the vec- 
tor components. We note that the momentum conser- 
vation equation and the spin evolution equation ([5]) 



still contains the explicit sum over the N states, and (as 
in classical fluid theory) it is necessary to impose fur- 
ther statistical relations in order to close the system [38[ . 
The preceding analysis applies equally well to electrons 
as holes or similar condensations. We will now include 
the ion species, which, due to the smaller charge-to-mass 
ratio, are described by the classical equations of motion. 

The coupling between the quantum plasma species is 
mediated by the electromagnetic field. By definition we 
let Bt t include spin sources, i.e. B to t = B + B sp , such 
that Ampere's law in terms of B to t reads V x Btot = 
A*o(j + jsp) + c~ 2 dtE, including the magnetization spin 
current j sp = V x(2nfiBS /%) 39]. We obtain consistency 
with the momentum conservation equation ([4]) by adding 
a term proportional to V X B sp to the Lorentz force, and 
subtracting it from the quantum force. The above alter- 
ations are only reshuffling of terms. However, a difference 
do appear when closing the system using Faraday's law. 
By letting V x E — —d t B totl using £? to t instead of B, 
we indeed obtain a difference compared to the classical 
Maxwell's equations. It is the full electromagnetic fields, 
including spin sources, that should be used in Faraday's 
law. Thus, Faraday's law as presented here is therefore 
the correct one to use. This form also gives a Hermitian 
susceptibility tensor (see below), something which is not 
obtained otherwise. 

To demonstrate the usefulness of the spin fluid equa- 
tions, we investigate linear wave propagation in a mag- 
netized plasma. For comparison we first neglect all 
quantum effects. Linearizing and Fourier analyzing the 
equations of motion, and substituting the velocities into 
Maxwells equations, we obtain e ■ E = 0, where 
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and the standard susceptibility components are 
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Here, the sums are over the particle species, k = (k 2 + 
k 2 ) 1 / 2 , fcx is the perpendicular (to z) part of the wave- 
vector, the T-direction is parallel to z x k±, lo p is the 
plasma frequency {Lo pe for the electrons and oj p i for the 
ions), to c = qBo/m is the cyclotron frequency, q and m 
are the particle charge and mass, v\ is the square of the 
thermal velocity times the ratio of specific heats, c is 
the speed of light in vacuum, and u)% — u 2 (w 2 — k 2 vf) — 
uj 2 (uj 2 — k 2 v 2 ). For notational convenience, the subscripts 
denoting the various particle species have been left out. 

Next we determine the equilibrium spin configuration. 
For many plasmas paramagnetic theory applies. Thus, 
in an external magnetic field B = B z, the zero order 
magnetization Mso due to the spin can be written [IoT ] 



Mso = n oVBV\-j^jr 



(10) 



where K is Boltzmann's constant, T is the tempera- 
ture, and we have introduced the Langevin function 
i](x) = [coth(a;) — x -1 ]. Here we have assumed that 
the spin contribution to the total magnetic field is small, 
otherwise Bo — > Bq + Bso, where Bso = MoMb^o, in 
Eq. (fT0|) . In general, the spin-magnetization Ms and 
the spin-vector S are related by 



S = 



hM s 



(11) 



and thus the zero order spin vector becomes Sq = 
(Ti/2)ri(iiBBo/KT)z. From (JTTJ) we obtain the spin- 
current contribution j s = V x (An e eS /m e ). 

Generalizing ([8]) to include all terms from quantum ef- 
fects give extremely complicated expressions. However, 
for most plasmas, the parameter /ibBq/ KT is very small, 
the spins are essentially randomly orientes, and the spin 
quantum effects are negligible. On the other hand, for 
low-frequency wave motion in a highly magnetized (or 
low temperature) plasma, the spin effects can be appre- 
ciable. In this case the dominant contribution to the spin 
effects come from the component of the spin force par- 
allel to the magnetic field, Fq z = — (2fiBnoSo/K)d z Bi, 
where B\ denotes the magnetic field perturbation, to- 
gether with the part of the spin current in the T-direction 
(from the part proportional to Vn x So), an d we can drop 
all other components as well as quantum terms that are 
proportional to h 2 , provided eBo 3> Tik 2 . Keeping the 
above terms, including only the lowest order contribu- 
tions in uj/uj c i, the susceptibility tensor is modified to 
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where the spin-contribution is 
^bBq 
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As an example we consider the fast and slow magne- 
tosonic mode, which is now described by the dispersion 
relation e TT e zz + (e T2 +£ sp ) 2 = 0. In the standard regime 
{jj 2 %/^K)i « 1, w « k z v te , the dispersion relation becomes 



{u) A -k 2 (? A )(u> 
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KT j m e Vf e 
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where the ion-acoustic velocity is c s — {m e /mi) 1 ' 2 Vt e , 
the Alfven velocity is ca = (B^/iionami) 1 / 2 and, for 
simpliticy, we have assumed that the ion-temperature 
is smaller than the electron temperature and included 
only electron thermal effects. Noting that ^bBq/KT = 
huj ce /m e v 2 e , obviously the spin-effects are important if 



huj r 



> 1. 



(15) 



Thus, for laboratory magnetic fields, where at most 
B ~ 10 — 20 T, we need low temperature plasmas for 
spin effects to influence the fast and slow magnetosonic 
modes. However, in the vicinity of pulsars and magne- 
tars [33|], we have B > 10 s T. For such systems, spin 
plasma effects can be important even in a high tempera- 
ture plasma. The spin effect on the fast and slow mode is 
illustrated in Fig. 1. Furthermore, we point out that for 
modes with even lower phase velocities (which exist in for 
example dusty plasmas [16]), the relative importance of 
the spin susceptibility term is enhanced, and spin effects 
can be significant also under laboratory conditions. 

In conclusion, we have derived the multi-fluid equa- 
tions for spin | quantum plasmas, starting from the 
Pauli equation. In order to demonstrate the usefulness 
of our equations, we have analyzed the linear modes, 
and demonstrated that the low-frequency modes are sig- 
nificantly altered by the spin effects provided the con- 
dition (fTSf is fulfilled. In many classical plasmas spin 
effects can be neglected due to the random orienta- 
tions of the spin vector. We stress here, however, that 
our results show that the spin multi-fluid equations can 
have important applications to such different mediums 
as low-temperature solid state plasmas, as well as to 
the accretion discs surrounding pulsars and magnetars. 
Furthermore, we emphasize that the spin-contributions 
are typically more imp ortant than the usual quantum 
plasma corrections 18J, specifically when the inequality 
eB > hk 2 is fulfilled. 

The linearized results presented in this Letter will most 
likely find its experimental application in dusty plasmas, 
where the low phase velocity will make the relative im- 
portance of the spin contribution (|13p particularly signif- 
icant, enabling probing of the collective spin dynamics. 

Finally, we suggest that the full nonlinear system ([3])- 
([7|) will show interesting behavior close to the electrons 
cyclotron frequency, when the spin vector evolution be- 
comes resonant. Moreover, the importance of the pres- 
sure like spin terms for, e.g., astrophysical plasmas is a 
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FIG. 1: The two roots of the dispersion relation (|14[) plotted 
in (a) (fast mode) and (b) (slow mode). In case (a) the lower 
surface is without spin and the upper surface is with spin, 
while in (b) the lower surface is with spin and the upper sur- 
face is without spin. We note that the contribution from the 
spin term can be significant, in particular for large values of 
the wave numbers k z and k±. Here we have used c 2 s /c 2 A — 0.5, 
rjhui ce /m e v 2 e — 8, and normalized the frequency by the ion 
cyclotron frequency Lo ci and the wave numbers by uj ci /cA- 



further field of investigation. 
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